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Abstract 



In case of an Ua{^) symmetry restoration in a hot and dense hadron matter, the 
mass of the produced hadrons and mesons can significantly change, and their 
production cross-section can also be modified. According to recent indirect 
measurement of Ref. [T] the rj' meson can suffer an at least 200 MeV mass 
drop. The decrease of the rj' mass results in an enhanced production, thus 
we can search for this effect via the decays of 1]' even if these decays happen 
in vacuum, outside of the hot and dense hadronic matter, where the mass 
modification happens. These decays of r/' change the two-pion correlation 
functions, and the dilepton spectrum. 

In this M.Sc. Thesis I search for the signature of an rj' enhancement in 
the PHENIX dilepton [2j spectrum in ^/suN = 200 GeV Au+Au collisions, 
which has a significant enhancement in the low-mass region as well, e.g. in 
the (0.1 — 1.0) GeV invariant electron-positron mass range. This is also the 
region of the t]' meson's Dalitz-decay {f]' — i- e+e^7), so the t]' enhancement 
might be responsible for at least a part of the excess. Other mesons' (other) 
properties can also be changed in the hot, dense medium, but in this thesis I 
focus on the mass modification of the rj' meson. To explore the role oirj', I used 
EXODUS [3] simulations to generate different dilepton spectra, corresponding 
to different rj' properties. The conclusion here was that the excess can not be 
described with just the 7]' mass modification, but the agreement with data has 
been improved significantly. 

Another idea which might brings us closer to understand the excess is to 
examine the radial flow of the mesons, which was not included in the original 
PHENIX analysis [2]. Radial flow is important in the low range, where 
it describes the particles' spectra well, just in the region where the dilepton 
spectrum has the most contributions from. Thus examining the effect of the 
radial flow seems to be inevitable, as it might be responsible for certain parts 
of the excess. 

The results summarized here are work in progress, internal results obtained 
with the framework of the PHENIX Collaboration at RHIC. I have presented 
them to the Collaboration in 3 talks during 2011-2012. 
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List of common abbreviations and notations 



BNL Brookhaven National Laboratory 

RHIC Relativistic Heavy-Ion Collider 

BBC Beam Beam Counter 

ZDC Zero Degree Calorimeter 

DC Drift Chamber 

PC Pad Chamber (PCI, PC2, PC3 correspond to the different layers of PC) 

RICH Ring Imaging Cherenkov Counter 

EMCal Electromagnetic Calorimeters, stand for both PbGl, and PbSc detectors 

Pt transverse momentum, Pt = Px ~^ Py 

rriT^ transverse mass, = \/p\ + m? 

y^SNN center of mass energy 

(ut) the average transverse flow 

To freeze-out temperature of the hadron gas 

Tefj effective slope of the of hydro spectra, e.g. 6"'"'^/"^'='^, where Tefj = Tq + m(uT)^ 



Chapter 1 
Introduction 



The results presented here are reports on a work in progress that proceeds in 
the framework of the PHENIX Collaboration. Published PHENIX data are 
presented from Ref. [2j for Au+Au collisions. 

Note, that this M.Sc. Thesis deals with one of the most difficult, unresolved 
and unsolved puzzles of relativistic heavy-ion physics, namely the low-mass 
dilepton enhancement in 200 GeV Au+Au collisions at RHIC. To our best 
knowledge, no published theoretical model has yet been able to describe these 
measurements. 
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Chapter 2 

The experiment 



2.1 RHIC (Relativistic Heavy Ion Collider) 

The experimental scene, the Relativistic Heavy Ion Collider (RHIC) is located 
at Brookhaven, in the center of Long Island, near New York in the USA. This 
collider is an intersecting storage ring of 3.8 km circumference, accelerating 
particles with alternating electromagnetic field. It is able to collide protons, 
and heavy-ions (Au+Au, Pb+Pb, Cu-Cu, U+U), and even particles with dif- 
ferent mass and charge, e.g. d+Au and Cu+Au. The acceleration energy can 
also be varied in a broad range, from ^snn = 7.7 — 200 GeV for Au+Au, and 
up to 510 GeV for p+p collisions, which makes it the most versatile collider 
ever built. It sets another world record by being the highest energy collider 
for polarized proton-proton collisions, but the corresponding spin physics will 
not be discussed in the current manuscript. 

The acceleration process is completed in steps, before the beam reaches the 
main ring, it goes through pre-accelerators (the Linac, Booster and AGS). An 
overview of the experiment and the acceleration setup is indicated on Fig. \2.M . 

Heavy-ions, using gold ions as an example, first go through the Tandem Van 
de Graaf linear accelerator, which uses static electricity to accelerate atoms re- 
moving some of their electrons. These ions leave Tandem with approximately 
5% of the speed of light (equivalent with 1 MeV energy/nucleon). The ion- 
ization here results Q = +32 e electric charge by the end of this stage. Later 
these ions enter two smaller synchrotons, the Booster and the AGS (Alternat- 



^ Figures [2TTII2. 41 are courtesy of the PHENIX Experiment at Brookhaven National Lab- 
oratory's Relativistic Heavy Ion Collider. 
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CHAPTER 2. THE EXPERIMENT 




Figure 2.1: Schematic picture of the RHIC accelerator with its pre-accelerators 

ing Gradient Synchroton), two circular accelerators, which accelerate them up 
to 0.37 c (95 MeV/nucleon) and 0.997 c (8.86 GeV/nucleon) and ionize them to 
Q = +77 e and +79 e electric charge, respectively. So actually before entering 
the RHIC ring, the gold atoms loose all of their electrons. 

Protons start in a different linear accelerator, the Linac, but after that they 
follow the same path as the heavy-ions. They leave the Linac with 200 MeV, 
the BOOSTER with 2 GeV and the AGS with 23 GeV energy. RHIC is also 
the only spin-polarized proton accelerator ever built, so its recent world record 
of s/s = 510 GeV for polarized p+p collision outdated an earlier record (500 
GeV) that was also reached at RHIC. 

After leaving AGS, the beam is passed into the AGS-to-RHIC transfer line, 
where a switching magnet splits the ions and protons to a so-called "blue" 
and "yellow" beams. At RHIC, two accelerator rings are located inside the 
tunnel which inherit the name of the beams injected into them. The blue 
one circulates ions clockwise, while the yellow ring accelerates ions counter- 
clockwise. RHIC's two concentric rings are made up of 1,740 superconducting 
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magnets. Due to further acceleration by RHIC, protons can reach ^Snn = 510 
GeV, the heavy-ions gain typically ^Snn = 200 GeV energy/nucleon. To op- 
erate effectively, the particles (both protons and ions) are collected into highly 
concentrated "bunches", each containing billions of ions, and 55-110 bunches 
circulate in a ring at the same time. 

The blue and yellow ring intercept each other in 6 points, 4 of which a detector- 
complex has been built on, the PHOBOS, BRAHMS, PHENIX and STAR 
experiments. PHOBOS and BRAHMS were designed as special purpose ex- 
periments, while STAR and PHENIX were foreseen as major user facilities 
with strong upgrade programs. PHOBOS was a tricky setup of silicon detec- 
tors, measuring charged particles from the collision, BRAHMS (Broad RAnge 
Hadron Magnetic Spectometer) was looking for a new state of matter, the color 
glass condensate (a detailed description of this state of matter can be found in 
Ref. jl]). Both detector-complexes were designed to detect particles with high 
pseudorapidities (close to the beam direction)|§, but their data taking period 
has already been completed: PHOBOS has been closed at 2005, BRAHMS a 
year later. The remaining two detector-complexes are STAR (Solenoid Tracker 
at RHIC), which covers the whole solid angle, and is able to reconstruct the 
path of charged particles with high accuracy, and PHENIX, which will be de- 
tailed subsequently, given that PHENIX research is the subject of this thesis. 
All the four RHIC experiments published the so-called "RHIC White Papers" 
in 2005, announcing the discovery of a new state of matter. More information 
about this discovery can be found in Refs. [SHS]. 

In 2010 PHENIX published its direct photon measurement in 200 GeV Au+Au 
collision, which revealed, for the first time experimentally, that the initial 
temperature of the hot and dense fireball and its nearly perfectly flowing fluid 
is hot enough to melt hadrons to a plasma of quarks (and supposedly gluons 
too). This measurement was based on virtual photon decays to lepton pairs, 
and the effective initial temperature was found to be Tinit > 300 MeV. The 
Hagedorn spectrum of hadrons suggests that hadrons as we know them can 
not exist at temperatures > Tnagedom ~ 180 MeV. 



^One definition of pseudorapidity is: 77 — ^ In jpj^p' , so high pz is equivalent with 
high pseudorapidity. The second one is directly related to the angle between the particle's 
momentum and the beam direction (here the z-axis): 77 = — In [tan |] . Comparing the 
two definition it is clear, that particles with high pseudorapidity fly out with small 9, cor- 
responding to particles that propagate close to the beam direction. 
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2.2 PHENIX experiment 
2.2.1 Introduction of PHENIX 

The PHENIX (Pioneering High Energy Nuclear Interaction eXperiment) 
detector-complex was designed to detect rare events of charged particles, and 
it is specialized to measure particles coming directly from the collision, the 
so-called direct probes of the quark-gluon plasma. These particles are not 
participating in the strong interaction, thus they can pass through the dense 
and hot medium without final state interactions, hence they carry information 
about the initial state and time evolution of the quark-gluon plasma. In par- 
ticular, PHENIX specialises on photons, electrons, muons as well as hadrons. 

PHENIX is a huge international collaboration, it involves ~ 500 scientists 
from 13 countries. It consists of several subdetectors, modules controlling these 
detectors and electronics responsible for communication and data-taking. The 
total weight of all three of the PHENIX magnets and steel are 1,657 tons, 
the whole weight of PHENIX is around 3,000 tons. Its dimensions are also 
impressive, it is 12 meters wide and four stories high. A schematic drawing of 
the PHENIX experiment as of 2012 is illustrated on Fig. 12. 2[ while a detailed 
cross-section drawing of PHENIX, showing the Central Arm and the Muon 
Arm can be seen on Fig. 12.41 

By location, its detectors can be grouped into the two central arms, where 
most of them are placed transverse to the beam, and the two muon arms which 
focus on the measurement of muons, that are located closer to the RHIC rings 
(north and south arms). By functionality, the detectors broadly fall into 3 
categories: the tracking detectors, the calorimeters and the event character- 
ization detectors. PHENIX also includes three huge magnets that bend the 
trajectories of charged particles. 

Tracking detectors measure the tracks of charged particles. The charged 
particles' momenta can be determined from curvature of their tracks in the 
magnetic field, if the PHENIX magnets are switched on. If these magnets are 
switched off, multiplicity and pseudorapidity density measurements are possi- 
ble with tracking. At PHENIX the Drift Chamber (DC) and the Pad Chamber 
(PC) are the most relevant tracking-type detectors. 

Calorimeters detect the energy of the particles by absorbing them or their 
decay products, this is why calorimeter-type detectors are always placed on 
the edge of a detector-complex, they are the last station for a particle. At 
PHENIX the Electromagnetic Calorimeters (PbSc and PbGl) and the Zero 
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Degree Calorimeter (ZDC) are the most important calorimeters. The PHENIX 
EMCals stop electromagnetic showers, but they are partially transparent for 
hadrons. 
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Figure 2.2: A shematic drawing of the PHENIX experiment as of 2012 



Particle identification detectors are based on various principles. For exam- 
ple, the Time of Flight detector (TOF) measures particle's speed obtained 
from their time of flight, and the path length of their tracks. Ring Imag- 
ing Cherenkov detector identifies electrons using their induced characteristic 
Cherenkov radiation. EMCals can also be used for particle identification, by 
measuring the energy of the particle with a given momentum, that is known 
from the tracking detectors. In the forward region special (MuID) plates were 
installed to identify muons. 

Beyond particle identification, to determine the reaction plane and the 
impact parameter of the collision is also fundamental, and this is a part of what 
event characterization detectors are responsible for. They also function 
as triggers, selecting interesting events and controlling the data-taking. The 
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reaction plane is determined by the beam axis and the impact parametecl, as 
illustrated on Fig. 




Figure 2.3: The reaction plane 

A short list of the PHENIX detectors, ordered by their location is given as 
follows: 

Detectors of the Central Arm: the Drift Chamber (DC), the Pad Cham- 
ber (PC), the Ring Imaging Cerenkov (RICH), the Time Expansion Chamber 
(TEC), the Time-of-Flight detector on the East side and the West side (TOF- 
E, TOF-W), the Aerogel Cerenkov Counter (Aerogel), and the Electromagnetic 
Calorimeters (PbSc and PbGl). 

The Muon Arm detectors: the Muon Tracker (MuTr), the Muon Identifier 
(MuID) and the Muon Piston Calorimeter (MPC). 

The Event Characterization Detectors: the Beam-Beam Counter (BBC), 
Zero Degree Calorimeters both on the north and south side (ZDC north and 
ZDC south, which are physically located inside the RHIC tunnel and their po- 
sition is not to scale on Fig. 12. 4p . the Forward Calorimeters, the Multiplicity 
Vertex Detector (MVD) and the Reaction Plane Detector. 



■^The impact parameter connects the center of the two coniding particles at the time of 
their closest approach. 
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West Beam View East 




I8.5m = 60 ft 



Figure 2.4: The cross-section drawing of PHENIX from different views 
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2.2.2 Detectors of the PHENIX dilepton measurements 

In this section an overview of the PHENIX detectors used in the dilepton 
spectrum measurement is presented. These are the global detectors; the Beam- 
Beam Counter (BBC) and the Zero-Degree Calorimeter (ZDC), two tracking 
detectors of the Central Arm; the Drift Chamber (DC) and Pad Chamber (PC), 
the Ring Imaging Cherenkov Counter (RICH) for electron identification, and 
the Electromagnetic Calorimeters (EMCal) for energy measurement, which 
also provide additional electron identification. 

The BBC and ZDC together measure the collision centrality and the vertex 
position along the beam direction. The centrality corresponds to the impact 
parameter of the collision, which is fundamental to measure. The ZDC and 
BBC are also used as first level triggers, so if they have hits simultaneously, 
data taking from all the other detectors can begin, otherwise signal in the de- 
tectors is considered as noise and data are not stored. The two BBC detectors 
are located 1.44 m from the collision point. They count the number of charged 
particles in pseudorapidity range 3.1 < \ri\ < 3.9, and deter the start-time of 
the collision with 20 ps resolution, which results in ~ 2 cm precision for p+p, 
and ~ 0.6 cm accuracy for Au+Au, respectively. 

The ZDC's are 18 m far from the interaction point in both directions al- 
most along the z-axis with a 3.6 cm vertical difference. They are focused 
on measuring neutrons coming from the Coulomb dissociation of beam nuclei 
and the evaporation neutrons emitted by the spectators. To understand the 
role of evaporation neutrons, consider a non-central heavy-ion collision. The 
peripheral part is not involved in the collision, it breaks off the nucleus and 
continues its movement along the z-axis. The neutrons and protons are almost 
uniformly placed inside the nucleus, so this splitting does not change the neu- 
tron/proton ratio. But bigger nuclei prefer to have more neutrons than the 
smaller ones, so this newly formed nucleus is not stable and emits a neutron 
with transverse momentum low compared to the typical 100 GeV momentum 
of the beam. These are called evaporation neutrons and due to the large beam 
energy used at RHIC, they move towards the ZDC. To be able to separate 
these neutrons from the spectators, other particles moving towards the beam 
direction and specially from the beam, the ZDC's were placed just after the 
dipole (DX) magnets, which bend the beam back to RHIC's two separate rings 
(the "yellow" and "blue" ring), thus keep charged particles off the ZDC's. 

The centrality of a collision is determined from these two detectors; it is 
defined as the correlation between the BBC and ZDC hits as seen on the top 
of Fig. 12.51 (Fig. from Ref. [9]). It shows the centrality classes determined 
of Run 2 PHENIX measurement |9]. The bottom part of Fig. 12.51 indicates 
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the impact parameters obtained from a HIJING [lOj simulation which has 
been filtered through the PHENIX specific PISA [llj software. Obviously, the 
centrality classes correlate with the impact parameters, but note that these are 
two different definitions. Note also that the two dimensional Qbbc/Qbbc ~ 
EzDc/ E^^Q plots allow for a more precise centrality determination than any of 
them alone: ZDC measures small number of neutrons both for most central and 
for most peripheral collisions, while BBC alone can not resolve the centralities 
below 40 - 50 %. 
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The Drift Chambers (DCs) and Pad Chambers (PCs) are in the Central 
Arm and are used to measure the transverse momenta of the particles from 
their bending curvature in magnetic field. The DCs are in the radial region 
from 2.02 m to 2.46 m. PC has 3 layers, the first (PCI) is between the DC 
and the RICH, PC2 is behind RICH (but only in the west arm) and PCS 
is in front of the EMCal. PCI and PCS are present in both arms. The 
magnetic field is produced by the Central Magnet which consists of two pairs 
of concentric coils able to run separately. It is an axial field magnet, and if 
both coils are running to the same direction, the momentum resolution for a 
single particle is better than 1% in the (0.2 — 1) GeV/c range. DC is applied to 
measure momenta in the plane perpendicular to the axis of collision, while PCI 
measures momentum along the collision axis, thus combining the two results 
in a precise momentum measurement. The additional layers provide pattern 
recognition and a possibility for systematic checks. 

The Time of Flight detector (TOP) measures the particles' speed from their 
time of flight in the following way. It has multiple tubes with length /, and 
every tube has two Photomultiplier Tubes (PMT's) at its both ends. If it has 
a hit at position x (measured from PMTl), the time required for the generated 
photons to reach PMTl is Ti = Tq + |, while the time they reach PMT2 is 
T2 = Tq + where v is the speed of photons moving in the tube, and Tq is 
the time it takes for a particle to emerge from a collision and to hit the tube. 
The starting time, namely the time when the collision happens is measured by 
the BBC's. Combining the two equations for Tq, and to take the arithmetic 
mean of them one gets the time of flight, tTOF = 



From the relationship txcF = ^ = ~ — p (-^ being the path length of a 

particle from the collision vertex to the TOF hit) one gets 



2 

m = — 



(2.2) 



where the momentum p is measured by the DC. Plotting the momentum as 
a function of m^, and multiplying it with the charge to separate matter from 
antimatter, one gets Fig. 
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2.2, 



Mass Squared [(GeV/c ) 1 



Figure 2.6: Particle Identification at PHENIX. The lines are the PHENIX 
specific cuts inside of which the signal is considered a pion, kaon and proton, 
respectively. Figure from Ref. |12| . 



This was just one method from the many used at PHENIX to identify 
particles. Another way is to calculate their mass from the energy measured 
by the EMCal and the momentum measured by the DCs and the PC's. For 
single electron identification a separate detector, the Ring Imaging Cherenkov 
(RICH) has been built, which is located in the radial region 2.5 m < r < 4.1 
m of the Central Arm. 
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2.3 PHENIX dilepton spectra 

Dilepton pairs are one of the direct probes of the quark-gluon plasma and the 
hot, dense medium formed in ^/s^ = 200 GeV Au+Au colhsions, since they 
can escape from the medium without final state interactions. By analysing 
the dilepton spectrum, one can probe the time evolution and dynamics of 
the collision, at least in principle. Dileptons can also be used to study the 
properties of low-mass vector mesons (p, u and (p) in the medium due to their 
relatively short lifetime compared to that of the medium (hadronic gas), and 
also to the pseudoscalar mesons, such as rj and rj'. Production of photons can 
be measured through their conversion to dileptons as well. 

The PHENIX Collaboration has measured the invariant mass spectrum 
(shown on Fig. 12. 8[ I2.10p . centrality (Fig. 12. lip and transverse momentum 
dependence (Fig. I2.12p of dilepton production in both Au+Au and p+p colli- 
sions [2] . This is a detailed investigation and we only highlight few parts which 
are the most important, and most relevant inputs of our study. PHENIX deter- 
mined the background dilepton spectrum from a p+p measurement. Different 
identified particles were measured, their spectra were fitted and parametrized. 
As for the particles' pt spectra, in the PHENIX analysis a modified Tsallis 
formula was used: 



E-^ = ATsams(e-("^-+^^^) + Pt/po)-" (2.3) 



With the formula above a simultaneous fit was elaborated to nearly all 
measured particles that decay to dileptons, as shown on Fig. 12.71 for p+p and 
on Fig. 12.91 for Au+Au collisions. 

Note, that proton and antiproton spectra has been also measured by TOF, 
but they are not included in these fits, since Fig. 12.91 contain only those that 
decay to dileptons. The proton+antiproton spectrum is important, however, 
since they have been measured in low transverse momentum range as well. 
As we shall demonstrate below, low particle decays dominate the min. bias 
200 GeV Au+Au dilepton spectrum, while the contributions of resonances with 
Pt > 2 GeV on Fig. 12.101 are nearly negligible. Thus it is very important to 
understand the low part of the meson spectra, and as seen on Fig. 12. 9[ 
already indicates some difficulties and incompatibilities with the PHENIX 
modified Tsallis function of Eq. 12.31 for Au+Au collisions, while in p+p colli- 
sions the spectra are well described by Eq. 12.31 as seen on Fig. 12.71 
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2.3.1 Dilepton in p+p 

As for p+p collisions, both the spectra and dilepton cocktail are well de- 
scribed as seen on Figs. \2.7\ 12.81 Note, that spectra also fit, which does 
not apply for Au+Au collisions, as detailed subsequently. 

Fig. 12.81 shows data and simulation results of Ref. [2J, elaborated with 
meson spectra of Eq. 12.31 This figure illustrates that for p+p collision 
the background is well understood. Given that in Au+Au collisions, however, 
both certain meson spectra and dilepton cocktail are not well described, this 
might confirm the presence of an r]' enhancement and/or a radial flow effect in 
Au+Au collisions, since in p+p no chiral symmetry restoration, nor radial flow 
effects have been observed so far. Hydro spectra with temperature gradient, 
however, might be at work. 




2 4 6 8 10 

(GeV/c) 



Figure 2.7: Different meson's px spectra of Eq. 12.31 for p+p collisions. Figure 
from Ref. 
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Figure 2.8: Dilepton cocktail of 200 GeV p+p collision, Figure from Ref. [2]. 
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2.3.2 Dilepton in Au+Au 

In the Au+Au, however, where both enhancement of Vj mesons and radial flow 
effect can play a significant role, the meson's spectra (as on Fig. I2.9p and the 
dilepton cocktail (Fig. I2.10|) are not completely understood yet. 

Note, that spectra, as on Fig. 12. 9^ are not well described by Eq. 12. 3^ 
possibly due to the lack of radial flow. An evidence for the radial flow effect 
in Au+Au collision and its absence in p+p collision is illustrated on Fig. 12.131 
of Ref. |13) . Note that the linear rise of the effective temperature with particle 
mass for non-central collisions was first derived in Ref. |24 |. 

As for the dilepton cocktail, a significant excess of measured dileptons as 
compared to the sum of known processes is seen in the (0.1 — 1) GeV region 
of Fig. 




Figure 2.9: Different meson's spectra of Eq. 12.31 for Au+Au collisions. 
Figure from Ref. 
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^ min. bias Au+Au \/s^ - 200 GeV 
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Figure 2.10: Dilepton cocktail of 200 GeV Au+Au collision, generated with 
meson px spectra of Eq. 12.31 Figure from Ref. |2]. 
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2.3.3 Centrality dependence 

The dilepton cocktail in different centrality classes has also been published 
in Ref. [2j, as seen on Fig. 12.111 This plot indicates the possible presence 
of hydro behaviour, as the excess tends to disappear for peripheral collisions, 
in accordance with the hydro-picture, where radial flow decreases with the 
centrality. Thus the excess can be a sum of (at least) two effects, one coming 
from a spectra which lacks the radial flow term, and the enhanced production 
of the Vj mesons, due to the in- medium mass modification effect. 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 

mee (GeV/c^) 

Figure 2.11: Dilepton cocktail of 200 GeV Au+Au collision, generated with 
meson spectra of Eq. 12. 3^ in different centrality classes. Figure from Ref. [2]. 
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2.3.4 Transverse momentum dependence 

Fig. 12.121 suggests that hydrodynamics-motivated spectra might result in a 
better agreement with data, since the input spectrum was fitted to the high 
Pt region, where the cocktail has negligible contributions from. The soft com- 
ponent of Eg. 12.31 lacks important physics ingredient, that we identify with the 
radial flow, thus it lacks to decribe the low-px part of meson spectra. If one 
looks at the meson peak of Fig. 12.121 at 1020 MeV, the low pt inaccuracy 
of Eq. 12.31 can be clearly observed, as the (f) meson peak has a different px 
dependence in the model/cocktail calculation as compared to the data. Also 
note that the low-mass dilepton enhancement is predominant in the soft, 
GeV < Pt < 1-5 GeV region, where the hydrodynamical picture is relevant. 
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Figure 2.12: Dilepton cocktail of 200 GeV Au+Au collision, generated with 
meson spectra of Eq. 12.31 in different p^ slices. Figure from Ref. [2]. 
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Figure 2.13: Inverse slope parameter (Tefj = Tq + m{uT)'^) obtained both for 
p+p and Au+Au collision. The fact, that for p+p it is almost constant if 
plotted as a function of mass indicates that radial flow effects (e.g. the average 
transverse flow, (ut)^) is present in Au+Au collisions but it is absent in p+p 
collisions. Figure from Ref. |13] . 
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Chapter 3 

Theoretical Introduction 



3.1 Chiral symmetry of the 3 quark model 

For setting the theoretical background, QCD's approximate chiral symmetry 
and references to chiral symmetry restoration models are outlined here. 

The term "chiral" comes from the Greek word for "hand", indicating that 
it is like a mirror symmetry between one's left and right hand. In particle 
physics, however, the chiral symmetry is defined as the invariance of the theory; 
the independent transformation of the left-handed and right-handed parts of 
Dirac fields in the QCD Lagrangian. Left-handed fermions correspond to a 
state when the particle's spin vector points to the opposite direction than its 
momentum. Right-handed state is when the momentum and spin vectors point 
to the same direction. 

To understand this phenomena in details, consider the QCD Lagrangian of 
the fermion Dirac-fields as an example: 

flavour flavour 

^f°= E W^^''D,)ql + qU^YD,)q],]~ ^ mHqlq], + q],ql) , (3.1) 

1=1 i=l 

where 7'^'s are the Dirac matrices, qi and q^ correspond to the left- and right- 
handed fermions projected from the fermion fields using = 0.5(1 — 75)5'*, 
and q)^ = 0.5(1 -|- 75)5* transformation. The q^ notation is the Dirac adjoint, 
e.g. Ql = q[i°- 

From Eq. 13. H it is easy to see, that the quark's masses explicitly break the 
chiral symmetry, e.g the mass terms are responsible for entangling the right- 
and left-handed parts of the Dirac fields, in accordance with the expectation 
from special relativity. According to that, massive particles can never reach 
the speed of light, thus if one considers a coordinate system moving any faster 
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than the velocity of a (e.g.) left-handed quark, viewing from that system, 
the momentum of that quark would point to the opposite direction, while its 
spin would be unchanged under this coordinate transformation. Thus, massive 
left-handed fermions can become right-handed after such a Lorentz boost. The 
quarks have to be massless to avoid the existence of such coordinate transfor- 
mations. So in the limit of massless quarks, chiral symmetry becomes an exact 
symmetry of QCD. So let's restrict ourselves to the massless quark approxima- 
tion {rriq = 0), and consider the SU{3) flavour symmetry, which contains only 
the u, d and s quarks and their antiquarks. In this case the transformation of 
the quark flavours will become unitary transformation, e.g. q\ — )■ U^l q'j^ and 
g)j — )■ f/]j g;^, where Ul-, Ur e SU{n), thus the SU{3) flavour symmetry of the 
quark model falls apart to a direct product of the left-handed and right-handed 
parts; SU{3) — )■ SU{3)l <S> ^^{'^)r- Thus, in the — > limit QCD becomes 
chirally symmetric. 

The chiral symmetry transformation has a different formulation, written 
in terms of the axial- vector symmetry. To do that, a new approach has to 
be introduced. Let Ul and Ur be the left-handed and right-handed unitary 
transformations defined as: 

f/^ = e^ei*", f/^ = e^®«*\ (3.2) 

Thus the quark Dirac fields transform as follows: 

The leading order, infinitesimal transformation reads as 

Introducing vector Qy and axial transformations, such as 



namely the vector part of this transformation which treats the left-handed 
and right-handed fermions equally, and the cixial part which treats them 
differently. Thus, the flavour transformation of a quark's Dirac field becomes 
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which means that QCD develops an approximate 

U{3)v ® U{3)A^SUi3)v ® SU{3)a ® Uv{l) ® Ua{1) (3.7) 

symmetry in the chiral hmit (e.g. in hmit of rUg — )■ 0), where f/y(l) is identi- 
fied by the barion charge and is conserved exactly, while the so-called Ua{^) 
symmetry is broken by the topological charge of QCD vacuum. 

This means that parity-partners have to have the same mass, which con- 
dition is generally not fulfilled, e.g. in the pseudoscalar octet m^r ~ 140 MeV, 
while mK,ri ~ 500 MeV. The difference between the quark masses would not 
result such deviation in the meson masses, only a few MeV difference would 
occur from that. It has been concluded, that the spontaneous breaking of the 
axial part ?7a(1) of the SU (3) symmetry is responsible for that mass deviation. 
In addition to that, the outstandingly large mass of the r]' meson (m^/ = 938 
MeV) complicates the picture even more. 

Further reading on the chiral symmetry restoration can be found in Ref . [23] , 
and also among its references. 
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3.2 The rj meson in a hot and dense medium 

In case of a chiral symmetry restoration in the hot and dense medium, the 
mesons can suffer a significant mass modification (and their other properties 
can also be changed). Here we focus on the %]' meson, and only on its mass 
modification effect. Its interaction cross-section is very small and its lifetime is 
very long, so it can escape from the hot hadronic medium without interacting 
with it, and will decay in the vacuum after regaining its original mass. Because 
of this, we do not have a direct observation channel for the chiral symmetry 
restoration through 77', but we have an indirect one by searching for enhanced 
production of Vj . So we distinguish two Vj mesons, and we note the properties 
of the in-medium 77' with an asterisk, while quantities without further notation 
correspond to the free 77' from now on. The Hagedorn formula estimates the 
enhanced production: 



where Tcond is the temperature of the condensate. This equation implies, that 
the lower the mass, the more 77' we have. This is the basic method applied in 
this Thesis, looking for signs of enhanced 77' production. 

Ref. [T] suggests a double exponential spectrum of r( in the transverse 
mas^, which predicts more particles to appear in the low region. This dou- 
ble exponential 77' spetrum will be used as an input in the current simulations, 
so an overview on the physical picture behind it will be given. 

So there are two components in the v( spectrum. Some of the v( have 
large enough transverse momenta so that they can escape from the hot and 
dense medium to the asymptotic vacuum by decreasing their momenta and 
increasing their mass, according to the conservation of energy: 



Those 77' mesons, however, that have low transverse momentunuin the chirally 
modified vacuum, cannot come on mass shell so they are trapped in that 
modified vacuum. They are released and come on-shell only when the modified 
vacuum state decays to the asymptotic vacuum. Both components realize an 
effective |0)^, + |?7% — )■ |0) + |r/') transition, but in a different way. 



^The transverse mass (tot) is defined from the transverse momentum, tot = \J to^ + . 



^Low transverse momentum means 77' mesons witli < \JmyY'^ — m*,^ as a result of 
conservation of energy (Eq. 13. 9p 




(3.8) 



* 2 I * 2 2 I 2 



(3.9) 
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The first component of the Vj spertum is the one coming from the conden- 
sate of which they were trapped into, to these a random transverse momentum 
was given: 



1 _ ''T 

/(Pt) = ^— e ^^V^. (3.10) 

27r m^/i? ^ 

Note, that most of the Vj mesons are emitted from the condensate at low px, 
and the transverse momentum generated in the Vj — )■ r)7r~^7r~ — > (vr+Tr'^Tr^) + 
7r+ + 7r~ chain is small, given the kinematics and the small difference of m^/ — 
5m^ ~ 250 MeV. 

The second, high part is a hydrodynamical spectrum: 

N{mT) = Cm!^e-"'T/^^«, (3.11) 

a = 1 — d/2, where d is the spatial dimension of the expansion (thus falls 
between -1/2 and 1). The effective temperature is Teg = Tq + m{u^)'^, where 
To is the freeze-out temperature and (ut)^ is the average transverse ffow, re- 
spectively. In this spectrum, rj' momenta are shifted according to Eq. 13.91 

The final rj' spectrum is the sum of Eq. 13.101 and Eq. 13.11^ which corre- 
sponds to a double exponential spectrum. In the current work the following 
parametrization was used: 

— = A ■ exp[— i?(mT — rriri')] + C ■ exp[— D(mT — "^r?')]- (3.12) 

The A, B, C and D parameters of the spectra and the enhancement 
(namely the r] and rj' weight factors) coming from this double exponential spec- 
trum were obtained in Ref. jl] for different resonance models (ALCOR |15] . 
Kaneta et al. [16], Letessier et al. |17] . UrQMD [18j and Stachel et al. |19] ) 
and are listed here in Table [3?T] for 77', and Table 13^2) for ry mesons. 

The input spectra for both rj' and rj mesons are shown on Fig. 13.11 The top 
left part shows the reconstructed (mx — to) spectrum of the rj' mesons using 
the resonance ratios of the Kaneta model [16j. The blue line corresponds to 
an 7]' spectrum without in-medium mass reduction, the red line shows the 
enhanced scenario. Top right part presents a comparison of reconstructed 
(mx — m) spectra of the rj' mesons for different models. The shaded (yellow) 
band represents the total error. Above — m'^ = 1 GeV, all models result 
in very similar values, corresponding to an approximate mx scaling. This 
figure indicates that the violation of this mx scaling is model dependent. The 
bottom left part of Fig. 13.11 shows a comparison of reconstructed (mx — to) 
spectra of the r] mesons from different models. The bottom right panel shows 
the absolute normalized spectra and compares it to PHENIX 200 GeV central 
Au+Au collision measurements 1221. 
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Model 


A 


B 


C 


D 




No enhancement 


0.82 


2.72 










ALCOR 


2.30 


2.98 


62.4 


23.5 


43.4 


Kaneta et al. 


2.21 


2.94 


32.4 


18.7 


25.6 


Letessier et al. 


2.91 


3.14 


80.1 


12.8 


67.6 


Stachel et al. 


2.85 


3.13 


80.0 


12.8 


67.6 


UrQMD 


2.75 


3.07 


52.5 


12.7 


45.0 



Table 3.1: rj' spectrum parameters of Eq. l3.12l from Ref. |T] for different models. 
The spectrum without enhancement was obtained from the Kaneta model 
using the m*, = m^/ substitution. 



Model 


A 


B 


C 


D 


/. 


No enhancement 


14.6 


3.38 










ALCOR 


14.6 


3.40 


97.0 


17.8 


5.25 


Kaneta et al. 


14.6 


3.38 


54.9 


16.2 


3.47 


Letessier et al. 


14.6 


3.38 


84.1 


16.9 


4.75 


Stachel et al. 


14.5 


3.38 


89.2 


17.0 


4.97 


UrQMD 


14.6 


3.41 


148 


17.9 


7.49 



Table 3.2: t] spectrum parameters of Eq. l3.12l from Ref. [T] for different models 
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Figure 3.1: r( and r/ enhancement as observed indirectly from tt^tt^ Bose- 
Einstein correlation measurements at RHIC, in 200 GeV Au+Au collisions, as 
in Refs. [Il|20] 
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Simulations of the 77 excess 



4.1 Introducing the simulation environment 

Our simulation tool, EXODUS is an event generator that has already been 
tuned to the PHENIX experiment [2]. It uses a fast Monte-Carlo method to 
simulate individual particles as well as their decays. It is capable of generating 
single particles, whole events, and most importantly, it can directly extract the 
dilepton spectrum in the specific PHENIX acceptance. Basically, this means a 
rapidity cut at \y\ < 0.35, and a Pt < 0.2 GeV/c cut in the transverse momenta 
of electrons and positrons, but further corrections have also been included. 

The program reads the particle information, e.g. the mass, width, charge and 
spin, from the "defined_particles.txt" file. For particle numbering it uses the 
standard Particle Data Group codes [H]. The decay definitions are in the "de- 
fined_decays.txt" file. One needs to specify the ID's of the parent and child 
particles, the number of child particles (NBody) and the branching ratio of 
the decay. Decays can be turned off by setting the option "Enabled" from 1 
to 0. Decay products ("children particles") can suffer further decays if "Chil- 
drenStable" is set to (and if that certain decay is defined). EXODUS can be 
run either interactively, or using an input file. 

In the input file, one has to choose between some predefined setups (e.g. 
"PHENIX single particles", "PHENIX: complete events" or the desired 
"PHENIX electron cocktail"). After some obvious steps, which include defin- 
ing the name of the output file (supports root and text format too) and the 
number of events to generate, the program requires the weight factor of vr^'s 
(which is basically the number of pions obtained from the integral of the 7r° 
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PT-distribution). After this, the weight factors have to be given for the other 
particles as well. In this Thesis minimum bias data was investigated, so the 
simulation has been done accordingly, by setting the number of binary N+N 
collisions to 257.8. At the end, the parameters of the pT-distribution (Eg. 12.31) 
have to be specified for each particle species, that decay to dileptons. PHENIX 
defaults to use a Tsallis distribution of Eq. 12.31 with the same parameters for 
all particles, except the normalization factor, Axsaiiis- The values of the pa- 
rameters are published in Ref. 
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4.2 Model predictions for the r]' and 77 enhance- 
ment 



The first estimation to understand the rj' meson's role in the dilepton spectrum 
was elaborated by using different input spectra from the PHENIX analysis |2] 
in EXODUS, for both the t] and t]' mesons, of which parameters have been 
obtained by an indirect measurement [Ij and listed here in Table 13.11 and 13.21 
In addition to the rescaling of r] and rj' with the obtained weight factors 
(according to the model predictions), additional t] mesons appear from the 
enhanced number of rj' mesons through the t]' — t- 71711] decay channel. Since 
none of the rj' ^ rj decays are calculated by EXODUS at this point, the rj 
weight factor has to be rescaled. The formula below naively estimates the rj 
mesons coming from the enhanced r]' mesons: 



where BR stands for branching ratio. 

So using the input spectra of Eq. l3.12l for rj' and rj, Eq. l2.3l for all the others 
an EXODUS simulation was elaborated for the different models. The simu- 
lation results for a/J/vw = 200 GeV Au+Au collisions with the x^/NDF and 
confidence level values in the region of the excess (0.15 GeV < Pt < 0.7 GeV) 
are shown on Figs. I4.m4.3[ Note, that this simulation has not been optimized 
yet to the best possible description of the dilepton spectrum, however, its fit 
quality for some of the models is already acceptable. 

This so-called first estimation obviously fails to explain completely the 
enhancement of the dilepton spectrum, but in some cases (e.g. for the Kaneta 
model) the results are quite promising. Given that, to our best knowledge, no 
other, similarly promising explanations up to the time of writing this thesis are 
publicly available, we follow up to this direction with more detailed simulations 
which might give a better qualitative-, perhaps even a quantitative description 
of the low-mass enhancement. 
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Original dilepton spectra 




ALCOR {r\ andri' enhancement) 




m„ (GeV/c") 

Figure 4.1: 

Top: EXODUS simulation result reproducing the dilepton spectra of the 
PHENIX analysis j2] without t]' enhancement, 
xVNDF (for 0.15 GeV < m^e < 0.7 GeV) = 4.6, CL = 0.011% 
Bottom: A promising first simulation including rj' and t] enhancement pre- 
dicted by the ALCOR [U] model, xVNDF (for 0.15 GeV < m^e < 0.7 GeV) 
= 2.25, CL = 1.8% 
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Kaneta {r\ andr|' enhancement) 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 




Figure 4.2: EXODUS simulation results including t]' and rj enhancement. The 
simulations are based on the results of [HEI], using resonance ratios from the 
Kaneta |16| (top) and the Letessier [T7] (bottom) models for multiplicities. 
Kaneta: xVNDF (for 0.15 GeV < TUee < 0.7 GeV) = 1.35, CL = 23.16%. This 
result is statistically acceptable. 

Letessier: x^/NDF is larger, than in the case without enhancement as it over- 
estimates the effect. 
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Stachel (r| andr|' enhancement) 




UrQMD (r| andr|' enhancement) 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 

m„ (GeV/c") 



Figure 4.3: Simulation result including i]' and i] enhancement predicted by 
the Stachel jTH] (top) and UrQMD [IB] (bottom) model. Obviously, these 
calculations are also promising, but they need to be fine-tuned to obtain a 
quantitative description of the low- mass dilepton enhancement. 
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4.3 Extending EXODUS: Simulation of decay 
chains 



So far only the direct dileptonic decays were included and calculated by EX- 
ODUS, other decays were taken into account via the fact, that inclusive px 
spectra were fitted which already contain the decays chain resulting the current 
meson. The normalization factors of the mesons are obtained from their px 
spectra and other independent sources (e.g. jet fragmentation) if data was not 
available. For further details, see Ref. ||2j. To explore the dilepton contribution 
of Vj in detail, however, specially when using different spectrum for that, 
the r( chain decays, e.g. Vj — ?■ riTT^n^ — t- (7r+7r°7r~)7r''7r'^ decay chains can play 
a significant role. But adding new decays slows down EXODUS significantly, 
and more importantly, one has to define the meson's weight factors at the be- 
ginning of the simulation (in the input file), which is not an effective way when 
elaborating a scan through different t]' weight factors, since its decay product's 
weights have to be corrected manually in each case. 

In this section a different simulation method is presented, e.g. the simu- 
lation is elaborated in each individual decay chain instead of generating the 
whole dilepton cocktail every time. This has the advantage that one is able 
to follow the contribution of individual particles, or simulate only the ones 
which parameters have been changed, and last but not at least this new ap- 
proach speeds up the simulation by a factor of 10 (as compared to its original 
speed). So in addition to the direct dileptonic decays already defined in EXO- 
DUS, other meson's decays were defined which result in a significant dilepton 
contribution, as listed in Table 14.11 

The consistency of this approach is checked via the rj' dilepton spectrum; it 
was first generated with all their decays defined at once (previous approach), 
then its individual decays were simulated and added up. The two results 
were cross-checked and are in accordance with each other, as seen on Fig. 14.41 
This indicates that using this method, generating individual decays instead of 
generating every decay of a given particle at once are two consistent methods 
and both are well suited to study rj' mass modification effects. 
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particle 


decay 


DIX 


iU 


^7331 


_j_ _ 

^211 ^-211 ^7221 
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n 
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0.29 
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0.2274 
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+ — 

e_ii 722 
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4 
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24 


7^111 


6-11 722 


n ni 7/1 
U.Ul ( 4 


o 
z 


<P333 


^7221722 


1.309 ■ 10-2 


115 




7r?ii722 


1.27- 10-3 


116 




6-11 ^ll 


2.954 ■ 10-^ 


117 




^22ieln en 


1.15 ■ 10"^ 


118 


J/^'443 


6-11 ^n 


5.94- 10-2 


119 


*444 


6-11 6ii 


9.7-10-6 


120 



Table 4.1: Table of the relevant decays contributing to the dilepton spectrum. 
The lower indices are the particles' ID's, according to the Particle Data Group 
convention [13], "BR" stands for branching ratio. "ID" is the decay identifica- 
tion number defined in EXODUS. Originally, only the mainly direct dileptonic 
decays were defined in EXODUS, while ID > 100 corresponds to the newly 
added decays and decay chains. 
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Reproducing thet)' spectra with individual decays 




0.2 0.4 0.6 0.8 1 1.2 

n\e {GeV/c") 



Figure 4.4: Reproducing the rj' spectra in two different ways: by simulating its 
relevant decays (ID 100-104 & ID 006) individually and adding them up (black 
line), and by defining the same decays at once (red line). The result shows that 
from individual decays the rj' spectrum can be reproduced faithfully. Thus this 
method will be applied to other mesons too, as its numerical implementation 
is much faster than that of the Monte Carlo simulation of each decay chains 
simultaneously. 



4.4 The m* — B ^ scan for r}' mesons 

So far only model predictions were used for the rj' spectrum, in this section we 
attempt to find the best rj' spectrum while scanning through the parameters 
of the 7]' spectra. 

The double exponential spectrum of Eq. 13.121 can be written in a different 
form (according to Ref. pj). The shape of the low-px part is governed by the 
temperature-like variable, -B~^, characterising the momentum distribution of 
the f]' mesons coming from the rj' condensate, the in-medium mass of rj' is 
responsible for the normalization (see Eq. 13. 8p . 

The and m* parameters define a certain rj' spectrum, so a scan through 
their physically reasonable values was elaborated with the two different ap- 
proaches. In the first one only the direct dileptonic decays are defined, and 
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this result will be contrasted to the case when all the r/' decays of Table 14.11 
are included. The second approach estimates the number of rj mesons coming 
from the enhanced number of r]' mesons better than the first one, the naive 
estimation of Eq. 14. ![ and changes the shape of the rj' dilepton spectrum as 
well. 

The maps of the results are on Fig. 14.111 Some dilepton spectra ex- 
amples for the best parameter obtained by this scan (e.g. for = 350 
MeV) for different mass values (m*,) are on Figs. I4.5M.7[ while Figs. I4.8M.10I 
show the dilepton cocktail for B^^ = 86.22 MeV (which was one of the best 
values of the previous analysis of Ref. [Ij), both with the same mass-drop 
values, respectively. These figures were simulated with the second approach, 
which involves every decays of Table H?T] for rj'. The original rj' contribution 
and the one coming from the enhanced number from its mass-drop are plotted 
separately. Note that because of the additional decays of rj', the shape of its 
dilepton spectrum has also been changed. 

The conclusion from these scans is that the rj' dileptonic contribution is not 
very sensitive to the slope of the low-px enhancement, the B~^ parameter 
as can be seen on Fig. 14.111 It is a conservative estimation, however, that 
this will play an important role when examining the cocktail in different 
slices. The important parameter is the in-medium mass of the rj', from the 
second approach it is expected to be around the rj mass. An other important 
learning from this section is that not only the direct dileptonic decays of rj' has 
to be included, as involving the resonance chain decays changes the dilepton 
spectrum, so that it gives a better agreement with PHENIX data. 

The right panel of Fig. 14.111 indicates that for any reasonable value of the 
slope of the rj' spectrum, B~^, at least 200 MeV mass drop in the medium is 
needed to describe the low-mass part of the 200 GeV Au+Au dilepton spec- 
trum, confirming the conclusion of Ref. [T], but based on a different channel of 
observation. This plot also indicates that it is important to take into account 
resonance decay chains of rj' properly. 



48 



4.4. THE M* - B-^ SCAN FOR i MESONS 





Figure 4.5: Examples from the result of the rj' input spectra scan for i? ^ = 350 
MeV, with mass values (from top to the bottom): 300 MeV, 450 MeV 
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I Electron Cocktail with excess 
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Figure 4.6: Examples from the result of the 77' input spectra scan for ^ = 350 
MeV, with mass values (from top to the bottom): 500 MeV, 600 MeV 
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Figure 4.7: Examples from the result of the rj' input spectra scan for = 350 
MeV, with mass values (from top to the bottom): 800 MeV and with the 
original mass, 958 MeV (in this case, no excess is expected) 
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I Electron Cocktail with excess 
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Figure 4.8: Examples from the result of the r]' input spectra scan for a more 
physical = 86.22 MeV, obtained from jl]. The mass values are (from top 
to the bottom): 300 MeV and 450 MeV. 
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4.4. THE M* - B-^ SCAN FOR i MESONS 





Figure 4.9: Examples from the result of the r]' input spectra scan for a more 
physical = 86.22 MeV, obtained from [1]. The mass values are (from top 
to the bottom): 500 MeV and 600 MeV. 
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I Electron Cocktail with excess 




I Electron Cocktail with excess 




Figure 4.10: Examples from the result of the v( input spectra scan for a more 
physical B^^ = 86.22 MeV, obtained from jT]. The mass values are (from top 
to the bottom): 800 MeV and with the original mass, 958 MeV (in this case, 
no excess is expected). 
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Figure 4.11: map for the two different approach. On the top, from sim- 
ulations with the 7]' meson's direct dileptonic decay-, on the bottom all the 
7]' decays of Table 14.1! are defined. On the vertical axis the parameter 
(which is responsible for the shape of 77' spectrum), on the horizontal axis the 
in- medium mass of rj' (m*,) is plotted, both in MeV units. The values are 
calculated for the region of the excess, e.g. from 0.12 GeV to 0.8 GeV, NDF = 
12. 
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Chapter 5 

The radial flow effect 



5.1 Why to use different spectra 



The px spectra of the PHENIX analysis ^ was fitted to the high-px region, 
as more data were available in that range. In the dilepton cocktail, however, 
unfortunately the low-pT part is dominant, for which a hydrodynamical spec- 
trum including radial fiow would be more a accurate estimation. But their 
formula (Eq. 12. 3p lacks the radial fiow term as shown on Fig. 15.21 To evaluate 
the importance of the contribution from low-px mesons, look at Fig. 15. where 
only dilepton pairs with (pT)ee < 2 GeV has contributed, hence, the original 
dilepton cocktail was reproduced only with them. Compare it with Fig. I2.10[ 
when the whole Pt range of the resonances up to 10 GeV was utilized. One 
can see that the addition of the 2 GeV < pi < 10 GeV region does not lead to 
extra enhancement. 

Another signature of the possible radial fiow effect is indicated by Figs. 12.111 
and 12.121 Fig. 12.111 shows the PHENIX cocktail in different centrality slices. 
The excess seen in the low-mass region tends to disappear for peripheral colli- 
sions, according to the hydro picture, where the (mt)^ parameter which governs 
the Pt distribution also decreases with increasing centrality. Fig. 12.121 shows 
the cocktail in different slices, and the introduction of a hydro spectrum 
is even more motivating here, since it shows that the excess is disappearing 
for higher transverse momenta. Thus the low-mass dilepton enhancement par- 
tially may be due to the lack of radial flow effect in the original background 
cocktail evaluation of Ref. |2j. A different input spectrum, which estimates 
better the low-px part of the resonances that decay to dileptons, may be an 
important tool in the upcoming systematic investigations. 
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Electron Cocktail p < 2 GeV 




m,, (GeV/c') 

Figure 5.1: To demonstrate the importance of the low-px part of meson spectra, 
only dileptons with (pT)ee < 2 GeV are plotted. They reproduce the dilepton, 
as this low-px part dominates the min. bias Au+Au dilepton cocktail. The 
simulation is new, while the PHENIX data points are from Ref. [2]. 



In the next section a more detailed examination of the PHENIX px spec- 
trum will be presented, and a dilepton cocktail generated with a hydro- motivated 
distribution will be shown, of which parameters have been obtained from a si- 
multaneous fit for TT^, K^, p and p. Due to the work-in-progress nature of the 
current results, in the next section we focus entirely on the effects of radial 
flow and the effects coming from resonance chain decays will not be considered 
in this section. For a final analysis, both radial flow and r/' chain decay effects 
have to be considered simultaneously, but these studies go well beyond the 
scope of this M.Sc. Thesis. 
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Spectrum of 



Spectrum of n* 




Spectrum of proton 




Spectrum of tj 





Spectrum of antiproton 




Spectrum of 




Figure 5.2: Comparison of the PHENIX spectra [2] (Eq. 12. 3p for different 
particles. The low-px approximation is the leading order of the exponential 
part of that spectra. It obviously lacks the radial flow effects, the mass-scaling 
inverse slope of the spectra, namely T^s = Tq + m{u^)'^ would imply that 
the more massive the particle is, the larger the tail in its exponential spectra 
(indicated by the dashed line), which is obviously not fulfilled here. 
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5.2. SIMULTANEOUS HYDRO FIT TO IDENTIFIED PARTICLE 

SPECTRA 



5.2 Simultaneous hydro fit to identified particle 
spectra 

As previous examinations suggested, a hydro-motivated spectrum might 
give a better estimation for the low-px part of the dilepton cocktail. The 
simplified hydro spectrum which is detailed in this section is: 

E—— = A( — ) e To+miut)-^ . (5,1) 
dp'^ \ m J 

The formula above was simultaneously fitted to available vr^, K^, p and p 
data, similarly to an earlier analysis of Ref. [,12], e.g. the fit ranges are (0.2—1.0) 
GeV/c^ for vr^ and (0.1 — 1.0) for K^, p and p in mx — m. The simultaneous 
fit to vr^, K^, p and p data is shown on Fig. 15. 3^ as a function of px- 




E I I I I t 3 1 I I I t 3 3 I I t I i I I I C ! 3 I I I I I I 

02 OA 0.6 0.8 1 1.2 1.4 1.6 1.8 



Figure 5.3: Simultaneous fit of Eq. 15.11 for vr^, K^, p and p data with the 
constraint from Ref. [12], e.g. the fit ranges are (0.2 — 1.0) GeV/c^ for tt^ and 
(0.1 — 1.0) for K^, p and p in — m. 



The best fit parameters are: {u^f = 0.21 ± 0.01, and Tq = (0.178 ± 0.004) 
GeV, while the a parameter was consistent with within errors, hence it was 
fixed to 0. Hopefully, this spectrum, when used as an input for generating of 
dilepton decays, and when extrapolated for heavy particles like r^', u and cf) 
will produce better agreement with PHENIX data. 
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5.3 Dilepton cocktail simulated with hydro spec- 
tra 



EXODUS simulations using Eq. 15.11 hydro spectrum as an input were elabo- 
rated, and are shown on Fig. 15. 4[ while Figs. l5.5l and r5.6l show its slices up to 
2 GeV. If Pt is beyond the hydrodynamical domain, and Eq. 15.11 is not appro- 
priate any more, the power-law tail QCD particle spectrum starts to take over 
then. Hence the p^ > 2 GeV region has been left out from this investigations, 
at present. 

The meson weight factors were obtained by fitting meson contributions as 
a linear combination to PHENIX data. It would have been better to obtain 
these weight factors from data, since the weight factors should come from the 
integral of their p^ spectra, but for all mesons contributing the dilepton spec- 
trum (p, u, (f), 7], rj' and even 7r°) no (or not enough) data were available in 
the pt < 2 GeV region, where the hydro formula of Eq. 15.11 is valid, so the 
normalization parameter {A of Eq. 15. ip could not be determined for them. The 
meson / pion yields coming from this fit are summarized in Table 15.11 Compar- 
ing 7]' /ir^ ratios, for current simulations a tenfold enhancement can be obtained 
if compared with the PHENIX analysis of Ref. p]. Using mesons as a basis, 
e.g. calculating rj' /(p ratios again for both cases results in even bigger enhance- 
ment. The physical picture behind this is that in the first case we divided the 
enhanced number of 7]'-s with an enhanced number of pions as well, as the rj' 
meson has a significant branching ratio to pions. The enhancement obtained 
from the comparison of ri'/(j) ratios is 38.1 for the rj', which is a same order of 
magnitude as of that reported in Ref. [T], based on indirect measurement of 
the 1]' spectrum. 

If one cross-checks the resonance ratios of Table 15.11 with Eq. 14. H will not 
get any contradiction, in fact, from that equation (using /^/ = 38.13) = 5.5 
can be obtained, while from the data (comparing r]/(f) ratios of Table EH] and 
Ref. |2]) an rj enhancement of 9.87 is resulted. The enhancements calculated 
from different methods are within a range of 1 magnitude. Possible additional 
effects, e.g. rj production cross-section might also be increased and of course, 
chain decays need to be more carefully simulated. 

Note, that the 1.5 GeV < pt < 2 GeV slice of Fig. 15.61 is fully explained. 
More work is needed to be done to combine t]' decay chains of Chapter 4 and 
radial flow effects of Chapter 5. 
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5.3. DILEPTON COCKTAIL SIMULATED WITH HYDRO SPECTRA 



meson 


meson/vr^ 


meson / 7r° of Ref . [2] 


V 


3.48 xlO-i 


L12 xlO-i 


P 


3.70 xlO-s 


8.98 xlO-2 


u 


L96 xl0~2 


L03 xlO-i 





6.74 xl0~3 


2.14 xlO-2 


v' 


2.57 xlO-i 


2.15 xlO-2 



Table 5.1: The meson/pion yield ratios obtained by fitting the different me- 
son contributions as linear combination. Note, that this method predicts a 
significant, tenfold enhancement of t]' mesons as compared to vr^-s. 



GeV < < 5 GeV 



T] -> e e y 
co-> e+e" & 71° e+e' 
e"'e' 

e^e^Y & e*e"e*e" 
71°^ e^e'Y 
(p^ e*e" & Ti e^'e" 
cc^ e*e (PYTHIA) 
sum with line 
200 GeV AuAu 




Figure 5.4: Dilepton spectrum simulated with the hydro distribution of Eq. 
15.11 The meson yields were obtained by fitting the meson contributions as 
a linear combination to PHENIX data. Note the significant increase of rj' 
production, as also shown in Table 13. Ij 13.21 and 15.11 As compared to neutral 
pions, a tenfold enhancement of rj' mesons appear as compared to the one a£ 
the PHENIX analysis [2]. 
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5.3. DILEPTON COCKTAIL SIMULATED WITH HYDRO SPECTRA 





Figure 5.6: Radial flow effects on the low-mass dilepton spectrum in different 
Pt slices 
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Chapter 6 
Summary 



The experimentally discovered low-mass dilepton enhancement is one of the 
most difficult puzzles of relativistic Au+Au collisions at 200 GeV. To date, 
this effect has not yet been described successfully. 

In this Thesis, I investigated if this dilepton excess in the low-mass region 
can be described by two independent effects; the enhanced production of the 
1]' mesons coming from the restoration of the chiral symmetry in the hot, dense 
medium, and the effect of the radial flow which has not been included in the 
PHENIX analysis (Ref. ^2]). 

As for including the radial flow term, ffist it was observed on Fig. 15. 
that only dileptons with < 2 GeV result in a significant contribution to 
the dilepton cocktail, thus finding a good description of the low-px part of the 
resonances' spectra seems to be inevitable. It was also observed on Fig. 15.21 
that the input spectra used in Ref. [2j (Eq. 12. 3p lacks the radial flow term, thus 
it estimates the (most important) low-px part of meson's spectra inaccurately. 
Hence, a simplified hydrodynamically motivated spectrum (Eq. 15. ip was in- 
troduced and fitted simultaneously to available vr^, ^ p and p data. The 
fit result for the freeze-out temperature (To) and radial flow term (ut)^ are 
in accordance with previous measurements and the fit resulted an acceptable 
confidence level of 8.6 %. The dilepton cocktail simulated with this spectra was 
then fitted to the dilepton cocktail of 200 GeV Au+Au collisions of Ref. [2]. 
The conclusion from this fit is that simulating resonances with a hydro Pt dis- 
tribution results in a good agreement with data, even for certain slices (1.0 
GeV — 2.0 GeV). Additional effects have to be considered as well to understand 
lower regions. 

As for the r]' excess, ffist dilepton cocktails with rj' spectra of Ref. [1] were 
simulated, as seen on Figs. 14.1114.31 These parameters were not fine-tuned 
to the current dilepton simulations, they resulted in a better agreement for 
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the low-mas excess for certain resonance models, however. Later, two scans 
through all physically reasonable r^' spectrum parameters {B~^ and m*) were 
elaborated, for simulations including only the direct dileptonic decays of rj' 
and with the extended version of the EXODUS event generator, when every 
significant decays of t]' were taken into account (the significant decays are listed 
in Table UT]). The map of the scan, for both methods is on Fig. 14.111 These 
studies gave no restriction for the slop^, but indicated an at least a 200 MeV 
mass-drop of the 7]', which is in accordance with the indirect measurement of 
Ref. [TJ. Note that the mass-drop of Ref. [Ij was elaborated in an independent 
observation channel, e.g. via tt^tt^ Bose-Einstein correlations. 

The examination of the radial flow effect resulted in an at least tenfold 
enhancement of the rj' meson (when obtained from r^'/vr^ ratio), such an en- 
hancement might also be a signature of the restoration of the chiral symmetry. 
When the rj' enhancement was obtained from comparing rj'/ ratios, however, 
an even higher, 40 times excess has been obtained (comparing to the PHENIX 
background from Ref. which is in the possible region of rj' enhancements 
obtained in Ref. [1]. 

We also found that proper simulations of rj' decay chains play an important 
role in obtaining the best values of the in-medium rj' mass, hence radial flow 
effects and r]' chain decays have to be combined, which is a subject of an 
upcoming investigation. 



^The quest for the slope parameter, B ^ might succeed, however, if dilepton spectra 
were fitted to different pT sHces. 
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